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Abstract 

Based on standard field-theoretic considerations, we develop an effective action approach for 
investigating quantum phase transitions in lattice Bose systems at arbitrary temperature. We 
begin by adding to the Hamiltonian of interest a symmetry breaking source term. Using time- 
dependent perturbation theory, we then expand the grand-canonical free energy as a double power 
series in both the tunneling and the source term. From here, an order parameter field is introduced 
in the standard way, and the underlying effective action is derived via a Legendre transformation. 
Determining the Ginzburg-Landau expansion to first order in the tunneling term, expressions for 
the Mott insulator-superfluid phase boundary, condensate density, average particle number, and 
compressibility are derived and analyzed in detail. Additionally, excitation spectra in the ordered 
phase are found by considering both longitudinal and transverse variations of the order parameter. 
Finally, these results are applied to the concrete case of the Bose-Hubbard Hamiltonian on a three 
dimensional cubic lattice, and compared with the corresponding results from mean-field theory. 
Although both approaches yield the same Mott insulator - superfluid phase boundary to first order 
in the tunneling, the predictions of our effective action theory turn out to be superior to the 
mean-field results deeper into the superfluid phase. 
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I. INTRODUCTION 

Recent developments in the field of dilute ultracold quantum gasses [1, 2, 3, 4, 5] have led 
to the experimental investigation of atoms in periodic potentials [6] . They are a fascinating 
new generation of many-particle quantum systems as they allow for the study of a variety of 
solid-state phenomena under perfectly controlled conditions [1, 6, 7, 8, 9, 10, 11, 12]. For in- 
stance, Bosonic lattice systems show a quantum phase transition for varying lattice depths. 
In deep lattices the tunneling between lattice sites is suppressed, and a Mott insulating 
state forms with a fixed number of Bosons residing on each lattice site. For shallow lattices, 
however, the dominance of inter-site tunneling allows for Bosons to coherently spread over 
the whole lattice, forming a superfiuid. The occurrence of such a quantum phase transi- 
tion between a Mott insulator and a superfiuid is observable, for instance, in time-of-fiight 
absorption pictures taken after switching off the lattice potential. They image momentum 
distributions integrated along one axis, and therefore by Heisenberg's uncertainty principle 
give information about the corresponding spacial distributions. Thus, the localization of 
atoms in the Mott phase results in diffuse absorption pictures, while the delocalized super- 
fluid phase gives rise to Bragg-like interference patterns. 

The theoretical analysis of this quantum phase transition is usually based on the Bose- 
Hubbard model Hamiltonian [13, 14, 15, 16], 



where cij and a] are Bosonic annihilation and creation operators, /i is the chemical potential, 
and signifies a sum over nearest neighbor sites i and j. Additionally, U parameter- 
izes the on-site interaction energy between two atoms at a given site, and t characterizes 
the kinetic energy, in this case given by the tunneling of an atom between two neighboring 
lattice sites. The quartic on-site coupling term, however, makes an exact diagonalization 
of (1) impossible. Thus, while Monte-Carlo simulations have proven fruitful for obtaining 
numerical results [17, 18, 19], analytic descriptions of Bosonic lattices near the quantum 
phase boundary have so far been typically limited to mean-field [13, 15] or strong-coupling 
approximations [20, 21]. Currently, the most precise analytic result for the whole Mott 
insulator-superfiuid phase diagram in a three dimensional cubic lattice at zero temperature 
is found in Ref. [22]. Therein, a Landau expansion for an effective potential with a spa- 
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tially and temporally global order parameter is derived. In this paper, we generalize the 
results of Ref. [22] by allowing for a spacially and temporally varying order parameter, thus 
determining a Ginzburg-Landau expansion for the effective action. This allows us to ob- 
tain an approximate analytic description of Bosonic lattice systems near the quantum phase 
boundary. 

To this end we proceed as follows. In Section II, we consider a very general type of 
Hamiltonian consisting of an arbitrary on-site interaction and an arbitrary tunneling term, 
of which the Bose- Hubbard Hamiltonian is a special case, and determine the grand-canonical 
free energy to first order in the tunneling term. This tunneling approximation is motivated 
by the fact that in three dimensions, the Mott insulator-superfiuid quantum phase transition 
is observed to occur for small values of t/U (note that the tunneling expansion is related to 
the random-walk expansion of Refs. [23, 24, 25]). Next, in Section III we introduce an order 
parameter field and derive a Ginzburg-Landau expansion of the effective action, allowing for 
the computation of physical quantities near the phase boundary in both the Mott insulator 
and the superfiuid phase. Sections IV and V present predictions of our effective-action 
theory for both static homogeneous and spatio-temporally varying order parameter fields, 
including expressions for the particle density, the compressibility, the superfiuid density, 
and the excitation spectra. Finally, in Section VI, we specify our results to the Bose- 
Hubbard Hamiltonian, and compare them to the predictions of the standard mean-field 
theory. Although both approaches yield the same approximation for the location of the 
phase boundary, our effective action approach turns out to be superior to the mean-field 
theory for the following reasons. First, we demonstrate that the effective action approach 
leads to qualitatively better results deeper into the superfiuid phase. Secondly, in contrast to 
the mean-field approximation, the effective action approach can be systematically extended 
to higher orders in the tunneling parameter in order to quantitatively improve the results, 
as has already been demonstrated for the case of the effective potential in Ref. [22]. 



II. GRAND-CANONICAL FREE ENERGY 

We consider Bosons on a background lattice with lattice sites denoted by i. Suppose they 
are described by a Hamiltonian of the form 

H = Ho + Hi, (2) 



which depends on Bosonic creation and annihilation operators at and Sj obeying the standard 
commutation relations 



0, 



Sij. 



(3) 



We assume that Hq is a sum of local terms each diagonal in the occupation number basis, 
i.e. 

Ho = J2Ma%) (4) 

i 

SO its energy eigenvalues are given by 



E 



(5) 



As we will be working grand-canonically, we stipulate that the terms fi{<x\ai) include the 
usual —fialdi dependence on the chemical potential /i. Also, we make the simplifying as- 
sumption that Hi is only a two-Boson hopping term 



Hi 



(6) 



with tij symmetric in i and j and = 0. As with the Bose-Hubbard model (1), Hq in Eq. (4) 
describes the Bosonic on-site interaction, while Hi in Eq. (6) incorporates the tunneling of 
Bosons between lattice sites. Note, however, that Eq. (2) with Eqs. (4) and (6) covers a 
significantly more general scenario than the Bose-Hubbard model. The on-site interaction in 
the Bose-Hubbard model (1) is a two-Boson term with a global interaction strength, but in 
Eq. (4), however, we have allowed for the on-site interaction of any finite number of Bosons. 
In addition, we have allowed the Hamiltonian to vary between lattice sites. Thus, our model 
is also capable of describing on-site disorder, which may arise from a local chemical potential, 
or from a local interaction [4, 13, 26]. Furthermore, Eq. (6) contains not only the tunneling 
of Bosons between nearest neighbor sites as in Eq. (2), but also between arbitrarily distant 
sites. 

As we are ultimately interested in investigating quantum phase transitions, we follow 
general field-theoretic considerations and add source terms to the Hamiltonian (2) in order 
to explicitly break any global symmetries [27, 28] 



(7) 
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Since the source currents ji{T),j*{T) depend explicitly upon the imaginary time variable r, 
standard time-dependent perturbation theory may be used to find a perturbative expression 
for the grand-canonical free energy. To do this, we switch to the imaginary-time Dirac 
interaction picture [29], with operators given by 



(8) 



where we have set ^ = 1. In this representation, the Schrodinger initial value problem for 
the time-evolution operator takes the form 

d 



dr 



Ud{t,to) = -H[-^{T)UB{r,To) 



UD{ro,To) 



This is solved by the Dyson expansion 



(9) 
(10) 



oo 

Ud{t, To) = l + ^Ui)' [T, To) 
n=l 



U^^\r,To) 



i-iy 



dri / dr2 . . . / dr^T 

TO J Ta J TO 



HUn)HUr2)...HUTnj 



(11) 

(12) 



where T is the standard imaginary-time ordering operator. The grand-canonical partition 
function for the system is defined as 

Z = tr{fe-^''^-^M}, (13) 

which can be rewritten as 

Z = tr|e-^^«f/D(AO)}. (14) 

This gives the partition function 2 as a functional of the currents. For brevity, we shall - 
in cases where no confusion may arise - suppress the arguments of functionals. Thus, by 
substituting Eqs. (11) and (12) into Eq. (14), we obtain 



n=l 



dri / dr2 . . . / dr„ ( T 

n'- Jo Jo Jo 



where 



zr=tr{e-'^^n=nE' 



-/3/«(n) 



(15) 

(16) 

(17) 



i n=0 



is the partition function of the unperturbed system, and 



< • >r 



1 



tr 



(18) 



represents the thermal average with respect to the unperturbed Hamiltonian Hq. This can 
be expressed more compactly as 



(19) 



Inserting the explicit form of H'i-[){t) from Eqs. (7) and (8), we see that the expectation values 
appearing in Eq. (16) can be expanded in terms of Green's functions of the unperturbed 
system. Furthermore, since the grand-canonical free energy is given as a logarithm of the 
partition function 



-^logZ, 



(20) 



the Linked Cluster Theorem [30] tells us that JF can be expanded diagrammatically in terms 
of cumulants defined as 



^2n (h 5 ''"l ! • • • ! '^n 5 '^U ■ ■ ■ ) 
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j=j*=0 



with the generating functional 



cr[j,/]=iog 



z 



Z(o) 



dr 





(22) 



with only contributions from connected diagrams [31]. Note that this approach, rather than 
a decomposition of the Green's functions via Wick's theorem, must be used in our case as 
Hq is not necessarily quadratic in the creation and annihilation operators. Because Hq is 
local according to Eq. (4), the average in Eq. (22) factors into independent averages for each 
lattice site. It follows that CQ^\j,j*] is a sum of local quantities, and thus the cumulants 
C^{i[, t[; . . i'^, t^|^i, Ti, . . r„) vanish unless all site indices are equal. With this, we 



can write 



(^l5 '^I'l • • • i ^n' ''nKl' ''"li • • • j Tn) — iiC*2n (''"l; • • • ; '^nl'^'li ■ ■ ■ i '^n) 



(23) 



n,m 



SO that it only remains to determine the local quantities iCf^{T[, . . . , r^|ri, . . . , r^). Using 
the definitions (21) and (22), we find that 



T 



6 



G'(°H^,nK,r2), 



(24) 



where G^^\i,Ti\j,T2) = SijG^^\i,Ti\i,T2) is the imaginary-time Green's function of the un- 
perturbed system. Similarly, 



Note that local the quantity iC^ (ri, T2\t3, T4) is symmetric under both the exchanges ri ^ T2 
and r-^ ^ r^. 

Because each power of the tunneling parameter tij is associated with a creation operator 
and an annihilation operator, and each power of jiir) {j*{T)) is associated with one creation 
(annihilation) operator, we can construct the connected diagrams which contribute to JF 
according to the following rules [32]: 

1. Each vertex with n lines entering and n lines exiting corresponds to a 2r2-th order 
cumulant iC2n • 

2. Draw all topologically inequivalent connected diagrams. 

3. Label each vertex with a site index, and each line with an imaginary-time variable. 

4. Each internal line is associated with a factor of tij. 

5. Each incoming (outgoing) external line is associated with a factor of jiir) (j *(''") )• 

6. Multiply by the multiplicity and divide by the symmetry factor. 

7. Integrate over all internal time variables. 

Each diagram is then multiplied by the appropriate factors of ji(T), jiir), and tij, and all 
spacetime variables are integrated. Since Hq in Eq. (4) is diagonal in the occupation number 
basis and local, there can be no contributions from diagrams with one line. Thus, to first 




^C^'\n\n),C^''>{T2\n) -.Cr(ri|r4).Cr(r2|r3). 



(25) 
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order in the tunneling tij and fourth order in the currents jiir) we find 



^ =Fo |/^^^ /^^^ ^2°^*^^' n\h ■r2)jiiri)j*{T2) + Ti\j, r2)tijji{Ti)j*{T2) 

I rP rP rP 
+ T / / dra / drg / dr^ a\ '{i, n; i, T2\i, T3; i, n)ji{Ti)ji{T2)fi{T:i)j*{n) 

4 Jo Jo Jo Jo 
I rP rP rP rP . 

+ / / dra / drg / dr^ V t.^ ar(i, n; i, T2\j, rg; i, n)ji{Ti)ji{r2)j*{T3)j*{n) 

^ Jo Jo Jo Jo J 

+ a'i\i, Ti; j, T2\i, n; i, T4)j,{Ti)jj{T2)j*{T3)j*{n) I , (26) 
where 

1 1 f °° 

e-/3/,(n) k (27) 



'"^ ^ i Ln,=0 



is the grand-canonical free energy of the unperturbed system, and the respective coefficients 
a2n are given by the following diagrams and expressions: 

af\z,n\z,T2) = n - > - T2 = ^Cf\n\T2), (28) 

a^\^,n\J,T2) = n » U » » r2= f dT,Cf\n\r),CP {T\r2), (29) 



r2 Ts 




af\i,Ti;i,T2\i,T3;i,T4,) = = iCfVi. ^2^3, r4), (30) 

Ti r4 




a^l\i,Ti;i,T2\j,T:i;i,n) = =/" dr^cf ^n, rsir, r4)jC5°Vk3)- (31) 



In the next section, we will see how this expansion (26) of the grand-canonical free energy 
leads to the Ginzburg-Landau expansion of the effective action to first order in ty. 

At this point, it is also worth making some observations about the two-particle Green's 
function G{i,Ti\j,T2), defined in the standard way as 



G{t,Ti\j,T2) = (f al{Ti)aj{T2) 



SMn)6jUT2) Z \ SUn)6j*{T2) 



(32) 
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This quantity can also be expanded diagrammatically in terms of cumulants, provided we 
realize that the effect of the prefactor 1/2 in Eq. (32) is simply to cancel all disconnected 
diagrams [29], ensuring that the only diagrams that contribute are connected diagrams with 
two external lines [31, 32]. Thus, there is a natural correspondence between the Green's 
function and the coefficients Og"'' defined above: 

G{i,n I J, T2) = (2, n | j, t^) + G^^^ (z, n | j, rs) + . . . 

= 5ij4°^(i, ri|i, T2) + (^5 Tili, T2) + .... (33) 



III. EFFECTIVE ACTION 



Evaluation of the diagrams shown in Eqs. (29) and (31) involve integration over the 
time variable associated with the internal line. Thus, their evaluation can be simplified by 
transforming to Matsubara space, where these integrals amount to simple multiplication. 
We use the following convention for the forward and inverse Matsubara transformations 

g{u:^) = l= [dre'--^g{T), (34) 



VP. 

^ 00 

^^(^) = ^ E aM^-'^"-^^ (35) 



m=—oo 



VP 

with the Matsubara frequencies 

27rm , . 

'^m = — ^ , m G Z. (36) 

Since the unperturbed Hamiltonian (4) is time-translation invariant, it follows from Eqs. (21) 
and (22) that the cumulants - and thus the functions - must depend on time-differences 
only. In terms of Matsubara frequencies, this implies that one of the frequency variables is 
restricted by a delta function, i.e. we have 

af\i, uJmi\i, uJm.2) = a2\h i^mi)Sujmi,iOm2y (37) 
Using this frequency conservation, we find from Eq. (29) 

a'i\i,UJml\j,UJ^2) = 4°^(«,t^ml)4°^(j,t^ml)^c^„l,a;„.2) (39) 
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and correspondingly Eq. (31) implies 



(40) 

Thus, the first order corrections to the functions a2n can be expressed entirely in terms of 
a~2^ and the corresponding zeroth order terms. Using the expressions (28) and (30), along 
with the definitions (24) and (25), we find that these two coefficients are explicitly given by 



^ oo 



n=0 



n + l 



n 



fi{n + 1) - fi{n) - iujm fi{n) - fi{n - 1) - iujn 



(41) 



and 



(0) / . . I . 

04 [1',^mU'i,i^m2\i,i^m4: 



^ 00 



n=0 



n{n — 1) 



fi{n - 2) - fi{n - 1) + iuJmA 

g/5(/i(n)-/j{"-l)-«^mi) — 1 \ 



+ 



^l3{fi{n)-fi(n-l)+i(Umi-^^ml~'^m2)) _ ]_ 

iiUmi - ^m2) \fi{n) - fi{n - 1) + iiuJmA " t^ml " t^m2) 

fi{n) - fi{n-l) - iLUmi ) fi{n - 1) - fi{n - 2) - iujrn2 
fi{n) - fi{n - 2) - i{ujmi + uJm2) fi{n) - fi{n - 1) - iuJmi J _ 

^2 r 1 / g/3(/i(n)-/i(n-l)+i{Wm4-'^ml-'^m2)) _ ]_ 



+ 



X 



fi{n) - fi{n-l) - iujm2 

g/3(/i(n)-/i{n-l)-ju;„ji) _ ^ 

fi{n) - fi{n-l) - icdmi 
1 

fi{n - 1) - fiin) + iuJmA 

n{n + 1) 
fi{n) - fi{n-l) - iuJrn2 



i{uJmA - t^m.2) \fi{n) " fi{n - 1) + i(u;m4 - t^ml " t^m2) 



^Pifiin)- fi{n-l)-ii^mi) — 1 \ 

/i(n) - /j(n - 1) - iw^i/ 
1 



/i(n + 1) - fi{n - 1) - z(t^„i + ujm2) 



oPUi{n)-fi(n-l)+i(u}m4,~^ml~^m2)) _ ]_ 



3/3(/i(n)-/i(n+l)+i(^,„4) _ ]_ 



/i('^) - /i('^ - 1) + K^mi - ^ml - ^m2) fi{n) - fi{n + 1) + iuJmA J 
I / g/3(/,(rj)-/i(n+l)+'t<^m4) _ I 



fi{n + l) - fiin) - iumi y^'^"'^™^'^™* /.(n) - fi{n + l) + iumi 

af^ (i, UJml\i, ^m4)af^ {i, UJm2) 



(42) 



where we have introduced the notation {•}r^^y to denote a symmetrization in the variables 
X and y. Hence the expansion of the grand-canonical free energy (26) can be compactly 
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rewritten as 



+ Nijkl{uJrnl,(^m2,(^m3,^m4)ji{(^ml)jj{^m2)jl{(^m3)jlii^rni) , (43) 



ijkl ^ml,l^m2 



where we have introduced the abbreviations 



(0), 



and 



]\T / \ ^i^ml+l^m2,l^m3+l^m4 (0)/- • I- \f X X X 



(44) 



kl 



tik(^2\^i ^m3)Sij + Uj(^2'{ji ^m2)^ik 



(45) 



Use of the expansion given above is hmited by the fact that the currents ji{uJm) are unphysical 
quantities. Therefore we desire a thermodynamic potential in terms of physically relevant 
observables. To this end, we define an order parameter field ipilum) in the standard field- 
theoretic way [27, 28] as 



5jT(w. 



(46) 



To first order in the tunneling parameter tjj, we find that the order parameter field is given 
by 



X Npjki{^ml,UJm2,UJm3,i^m)jp{uJ,nl)jj{uJm2)jl{uJm3)- (47) 

pjk iOml,^m2,i^m3 

This finding motivates the performance of a Legendre transformation of JF to obtain the 
effective action which is a functional of the order parameter field: 

T[i)i{u;m),'iPt{^m)] = ~ ^ X X [^i(^rn)3l{^m) + i^* {^m) 3i{^m)]- (48) 

i LUrn 

The importance of the functional F is made clear with the following observation. The 
physical situation of interest is the case when the artificially introduced currents vanish, i.e. 
when we set iiitOm) = ilij^m) = 0. Since ip and j* are conjugate variables, we have that 

5T 



ji{uJm) = -P 



m I 



(49) 
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and thus this physical situation corresponds to 



6r 



6r 



1p=1pe: 



0. 



(50) 



?/)=l/lc 



This means that the equihbrium value of the square of the order parameter field is 
determined by the condition that the effective action F is stationary with respect to variations 
about it. Furthermore, we have from Eq. (48) that the effective action F, evaluated at the 
equilibrium order parameter field, is equal to the physical grand-canonical free energy: 



lim JF. 



(51) 



Now, a Ginzburg-Landau expansion of the effective action can be obtained. First, using the 
fact that to first order in 



5„ 



m2 



Eq. (47) can be inverted recursively to find ji{oJm) as a functional of the order parameter 
field, yielding 



6ij — a'^\i, ujmi)ti 



(52) 



^p(^ml) 



qjk I^m2,l^m3 



(53) 



where we have defined the abbreviation 



(54) 



Inserting this expression for jiioJm) into Eq. (48) together with the expansion (43), and 
keeping terms only up to first order in the tunneling tjj, we find 

|2 



F = Fr 



(55) 



(0) / . . I . . X 

^ A (0)/- N (0)/. N (0)/. N (0)/. r 

l^5'^mlja2 {^■y'^m2)0'2 ['''y^m3)C'2 l^! '^m4j 



Thus, after performing the Legendre transformation, it turns out that the tunneling param- 
eter tij appears up to first order only in terms which are quadratic in the order parameter 
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field. Furthermore, note that this result for the effective action is sufficiently general that it 
depends on only three quantities of the unperturbed system: the grand-canonical free energy 
(27) and the Matsubara transform of the zeroth-order coefficients (41) and (42). Finally, 
the condition for equilibrium (50) becomes 







(0)/- 



E 



,(0), 



2a 



(0), 



A{l^ml)A{^m2)^i{l^m3)- (56) 

Due to the complexity introduced by allowing the functions in the Hamiltonian (4) 
to be site-dependent, and the fact that many interesting physical scenarios can be modeled 
with a uniform on-site interaction, we restrict our attention in the rest of this paper to the 
homogeneous situation 

fiiajcLi) = f{aldi). (57) 

In this case, the cumulants are no longer on-site quantities, and we may thus drop the site 
indices in the coefficients a2n- In the next sections, we examine the physical implications of 
both a static and a dynamic order parameter field. 



IV. PHYSICAL QUANTITIES IN THE STATIC CASE 



Consider first an order parameter field that is constant in both time and space, i.e. of 
the form 

tpi{uJm) = '4'VpSum,0- (58) 

With this, the effective action (55) simplifies to the effective potential 

/3af (0,0|0,0) , ,,4 



af(0) 4 



af(0)' 



7 + ^0 



(59) 



where A^^ denotes the total number of lattice sites, and 7 = ^j^ tij- In the case where 

af^(0,0|0,0) < (60) 

we have, according to the standard Landau theory, a phase transition of second order with 
a phase boundary given by the set of system parameters satisfying 

AT, 







-7. 



(61) 
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As this is the case of most interest, we will assume that such a phase transition exists. We 
also find that Eq. (56) takes the simple form 

A^. . ,.2PN,af\0,0\0,0) 



= ij 



from which we see that in the ordered phase the equilibrium value of 
condensate density, is given by 



(62) 



eq 



and thus the 



(63) 



/5iV,af (0,0|0,0) 

Furthermore, due to Eq. (51), other physical quantities follow from evaluating derivatives 
of r at ipeq. For instance, the expectation value of the number of particles per lattice site 
{n) = in the ordered phase is given by 



1 dV 

in) = -TT- 



(64) 



'0=V'cq 



and correspondingly, the compressibility k, 



follows from 

Ofl 



1 d^T 



(65) 



In general, any thermodynamic quantity expressible as a function of derivatives of the grand- 
canonical free energy JF can be expressed as the same function of derivatives of F with respect 
to the same variables, evaluated a.t = ip^q. 

Relaxing the condition of spacial homogeneity of the order parameter, we are able to 
determine the superfiuid density of the system. The superfiuid density is defined as the 
effective fluid density that remains at rest when the entire system is moved at a constant 
velocity [33, 34]. As is well known in quantum mechanics, such a uniform velocity corre- 
sponds to imposing twisted boundary conditions. Equivalently, we introduce Peierls phase 
factors 

cti CLic'^-^ (66) 

in the original Hamiltonian (2). Here is related to the velocity of the system according 
to t? = (f)/m*L where m* is the effective particle mass, Xj are the lattice vectors, and L is 
the extent of the system in the direction of v. Equating the kinetic energy of the superfiuid 
with the free energy difference ^(0) — ^(0), we see that the superfiuid density p is given by 



lim — 

A^,|0|^ 



P 



^(0)-^(O) 



(67) 



14 



Examining the form of Hq and Hi in Eqs. (4) and (6), we see that the effect of introducing 
the phase factors in Eq. (66) is simply to redefine the tunneling parameter tij as 



Thus, using Eq. (51) we can express p in terms of the effective action as 



p = lim 



which, with the aid of Eq. (55) reduces to 



V' = '/'oq(0) 



r(o) 



»/'=V'cq{0) 



p = lim 



u 



V^cq(0) 



2 



e L - 



^cq(0) - ^cq(O) 



^eq(O) 

/3af^(0,0|0,0) 



^eq(0) - ^cq(O) 



(68) 



(69) 



(70) 



4^0) V^^^^-^ J 4(af(0))^ 

Thus, the superfiuid density is determined explicitly once a definite form of the tunneling 
parameter tij is specified. 



V. PHYSICAL QUANTITIES IN THE DYNAMIC CASE 



By allowing the order parameter to vary in imaginary time, we can also use the effective 
action to obtain an analytic form for the Matsubara Green's function. To do so, we note 
that the Legendre transformation (46), (48), (49) implies 



(3 



6^T 



ml ) 



(71) 

We recognize immediately from Eq. (32) that this is precisely the inverse of the Matsubara 
Green's function Q{i,uJmi\j,^^m2)- Next, we consider T expanded to arbitrary order in Uj, 



^ i Vo;™ 4 (^m) n=l j 



+ Y.Y.''^^^M[miAMr,M + ...], (72) 



where the expansion coefficients a^'^ are determined by methods like those described above. 
We then find from Eq. (71) that the Matsubara Green's function is given by 



(73) 



n=l 
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Recognizing that Suj^^^ujm2^ij / o,^\^mi) is simply the inverse of the unperturbed Matsubara 
Green's function, we see that the power series in tij in Eq. (73) gives a series expansion of 
the self-energy S: 

oo 
n=l 

(74) 

Thus, we conclude that our effective action gives an expansion for the Green's function in 
terms of the self-energy in powers of tij. In the non-ordered phase, this is the same as if 
we had computed the corrections to the unperturbed Green's function directly from our 
perturbative expansion of JF and performed a resummation [32]. Specifying Eq. (73) to our 
present first-order case, we hence find 



Si 



26, 



(a;^i,0|0,a;^2)(ar(0))M ^ 



(0), 



(^2^ {^ml ) 



+ 



4«)(o) 



(af (^^i))2af (0,0|0,0) 



tij. 



(75) 



Denoting by ttt, the Fourier transform of the tunneling parameter. 



tijC 



(76) 



we find that Eq. (75) can be rewritten in Fourier space as 



g{k',U,n2\k,UJml) 



1 



25. 



(a;„,i,0|0,u;^2)(af (0))2' 



+ 



t 



kk> 



(77) 



afM (4°)(..„i))2af (0,0|0,0) 

We see that the second term above is a contribution to the Green's function due completely 
to the existence of a non-vanishing order parameter. This correction can thus be exploited 
to improve analytical time-of-fiight calculations for Bosonic lattice systems in the superfluid 
phase [35]. 

Next, we can examine excitations of the system at zero temperature by looking for spatio- 
temporal variations of the order parameter field about ip^c^ = \/\i^\1q^^^° which preserve the 



equilibrium condition (50), where 6q is an arbitrary global phase. To this end, we first must 
specify to the case where our system is translationally invariant, such that 



k"kk' 



(78) 
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Next, we add to the equilibrium value of the order parameter field a small variation 
6'il){xi,uJm)- We then Taylor expand the effective action F about ip^q in terms of these 
variations. The first order term vanishes due to the equilibrium condition (50), leaving 



6^r 



6{6llJi{uJml))S{Sllj*{uJm2)) 



m2, 



+ ... 



(79) 



Demanding that in equilibrium the effective potential is stationary with respect to the 
variations dip in the standard way gives the equation of motion 

6iJi{Uml) = 0. (80) 



6^r 



Stp=0 



^ 5(5^i(a;„,i))5(5^*(^^2)) 
This equation can be satisfied in two distinct ways. The trivial solution 6ip{xi, Um) = 
corresponds to the static homogeneous equilibrium examined in the previous section. The 
second solution is given by 

= (81) 



5ip=0 



and describes the excitation spectrum of the system. In particular, by analytically contin- 
uing Eq. (81) to real frequencies and transforming to Fourier space, we are able to identify 
the dispersion relation of low-lying excitations as those curves uj{k) which make the equa- 
tion valid. The standard method of performing this analytic continuation is to find the 
equations of motion in imaginary-time and perform an inverse Wick rotation. Because of 
the complexity of the coefficient a^4\uJmi,uJm2\i^m3,i^mA), however, this is ill-suited to our 
present needs. Therefore, we note that our imaginary time evolution operator exp (^—Ht^ 
can be mapped to the real-time evolution operator exp (^—iHt^ by the formal substitution 
H iH. To maintain the reality of the grand-canonical free energy, we must also perform 
the substitution JF —iT. We thus find that in terms of real frequencies the effective 



action is given by 



i dui dLJ2 dujs 




YtijilJi{uj)ilj*{ 



(82) 



■ iafl {ui ) a^l {lj2 ) a^2R (^s) 4r (^4) 
where and obtained from Eqs. (41) and (42) respectively by the replacement 

fi{n) — > ifi[n). Thus, the real-time continuation of the condition (81) is given by 



^2Ft 
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0. 



(83) 



In general, the function u!{k) will have a positive and a negative frequency branch. Because 
we determined these curves by expanding the effective action Fr about a minimum, however, 
only the positive frequency branch of uj{k) are to be considered as physically relevant. 

Since the order parameter is complex, we examine separately variations of both the 
magnitude and of the phase. First, we consider excitations in the amplitude of the order 
parameter. To this end, we replace ip in Eq- (55) by ipcciVPSuj„,,o + Si^i{^^m), where 6ipi{uJm) is 
an arbitrary infinitesimal function of the lattice site i and Um., with fixed phase 6q. Carrying 
out the functional derivative in Eq. (83) and performing the continuation outlined above 
and transforming to Fourier space yields the equation 



2aiS(^A,0|0,^A)(4°R(0)) 







+ 



I + aiiO) 



(84) 



4!;(^a) ' (aS(o;A))^aS(0,0|0,0) 
This gives a constraint equation which can be solved for the dispersion relation of amplitude 
excitations ujx{k). By comparing Eq. (84) with the Matsubara Green's function (77), we 
notice that the dispersion relation ujA{k) coincides with the poles of the translationally 
invariant real-time Green's function. 

To treat the phase degree of freedom, we note first that adding a small time- varying phase 
to ipeci amounts to the transformation 



Expressing this in Matsubara space, we have 



eq 



(85) 



(86) 



Inserting the transformation (86) into the real-time effective action (82) and performing the 
derivative (83) yields the condition 







—I 



2aS(0)' 



+ 



+ 



4r(0) 



a^liug) 42(0) 42(0,010,0) 

+ 6(0, 0, Ug, -UJg) - 2b{ujg, 0, Uq, 0) - 2b{uJg, 0, 0, Ug) 

where we have defined 



26(0, 0,0,0) + 6(cue, -cue, 0,0) 



(0)^ I ^ 



4S(wi)4r (^2)4r (t^3)4R (^4) 



(0), 



(87) 



(88) 
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This determines the dispersion relation ujg{k) of the phase excitations. We note that since 
% = l/Ns, oJeiJy) = is a solution to the constraint (87) in accordance with Goldstone's 
theorem [27, 28]. 

Lastly, we investigate the phenomenon of second sound. As is well known, the observed 
elementary excitations of a superfluid are given by phonons. To obtain their corresponding 
dispersion relation u!s{k), we must examine the phase excitations in the presence of the am- 
plitude variations, i.e. ip [ipcq + ^'^i(^^m)] e*^'*^"^'"^. This has been considered, for example, 
in Refs. [36, 37], leading to the result 



uj,{k) = ^ujA{k)uJe{k). (89) 

VI. AN APPLICATION: THE BOSE-HUBBARD HAMILTONIAN 

Having developed the field-theoretic approach for the general Hamiltonian (2), (4), (6) 
in the previous sections, we are now in a position to apply it to the specific case of the 
Bose-Hubbard Hamiltonian on a 3-dimensional cubic lattice defined by Eq. (1). As is well 
known, this model exhibits a quantum phase transition between a Mott insulating phase 
and a superfluid phase [13, 14, 15, 16]. The Hamiltonian (1) has exactly the form assumed 
in Section II when the following identifications are made: 

f{n) = En = Ujn{n - 1) - /in (90) 



where cr^ denotes the lattice basis vectors with A; = 1, 2, 3. Additionally, we see that the 
quantity 7 introduced above simplifies to 

l = Y.Uj = QNst. (92) 

Thus, we can apply all previously derived formulas to extract physical information about 
the Bose-Hubbard system. In the Um — * limit, we find that a2°''(0) becomes 



The expression for a4'^''(0, 0|0, 0) must be calculated via a more careful limiting procedure. 
Making use of the limit 

lim = 6, (94) 
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we find that 



af (0, 0|0, 0) = g e-^^" - 1) 



+ n 



+ 



/5 



2/9 



(-E'n — En-lY{En — En-2) 

2(-E'n+i — 2i?„ + -E„_i) 



+ 1 



(-E'n+l — En){En — En^i) _ 

+ (n + l)(n + 2) 



n + 1 



/3 



.(-E'n+l — EnY{En+2 — En) _ 



[{En+i-Enf {En+i-Enf 



(95) 



A. Effective Action Predictions in the Static Case 



Before examining the physical imphcations of our effective action approach, we first intro- 
duce the standard mean-field treatment of the Bose-Hubbard Hamiltonian for comparison. 
The mean-field Hamiltonian is found by performing a Hartree-Fock expansion of the hopping 
term in the Hamiltonian (1) [13, 15]. Keeping in mind that the order parameter is defined 
according to ip = (Sj), this yields 

^MF = i^o - 6t ^ (^l)a\ + - IV^I') . (96) 

i 

The methods of Section II can be adapted to give an expansion of the mean-field free energy 
in powers of the order parameter, since 

^MF = H^ + H[\t=o + QtNs li^f , (97) 

when we make the formal identification jj(r) = —Gttp. Thus, an expansion of JF to zeroth 
order in t gives an expansion of the mean-field free energy J-'mf in powers of the order 
parameter, provided we recognize that the constant term in Eq. (97) contributes a term of 
order to J-'mf- With these considerations in mind, we find the explicit result 

^MF = F,-N, \^\' + far l^l') , (98) 

where the mean-field Landau coefficients a^^ and a^^ are given by 

af"" = af\0)iQty - Qt, (99) 
= af (0,0|0,0)(6t)^ (100) 
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FIG. 1: Plot of the critical value of the hopping parameter t versus the chemical potential /i, 
both scaled by the interaction energy U. The phase boundaries for two different temperatures are 
shown. The solid blue curve \sT/U = 0, and the dashed red curve is T/U = 0.1/kB- 

Thus, the mean-field result can also be expressed in terms of the same three quantities (27), 
(41), and (42) as our effective action approach. 

We now compare the predictions of our mean field theory with our effective action theory. 
First, we find that the mean-field phase boundary is given by the curve [13, 26, 38, 39] 



which turns out to be identical to the phase boundary found from Eq. (61). A plot of 
the phase boundary in Fig. 1 reveals that increasing thermal fluctuations destroy quantum 
coherence, as the superfluid phase shrinks with increasing temperature. Note that the main 
advantage of the field-theoretic method over the mean-field approach is in the fact that the 
phase boundary can be improved by carrying the expansion out to higher orders in t. The 
phase boundary to second tunneling order has already been calculated for T = in Ref. [22] 
and for T > in Ref. [32], and proves to be a considerable improvement over the mean- field 
result. 

We next look now at the condensate density {ipll^- From Eq. (63), we find 



.MF 



1 



(101) 





(102) 



while standard Landau-theory yields 




(103) 
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FIG. 2: Plot of the condensate density as a function of the hopping parameter t/U for both the 
mean-field theory (left) and the effective action theory (right), with fixed fi/U = 0.9. The solid 
blue curves are T/U = 0, and the dashed red curves are T/U = 0.1/kB- 

via the minimization of J-'mf- 

Turning to the superfluid density, we see from Eq. (68) that for the Bose-Hubbard model, 



where d^^ are the nearest neighbor lattice vectors in the a direction. Therefore, in the Bose- 
Hubbard model we have t = l/(2m*) [10]. Thus, from Eqs. (67) and (70) we find that for 
both the mean-field theory and the effective action theory to first order in t, 



Taking to be in a lattice direction, we see that the superfluid and condensate densities are 
equal at this level of approximation. 

A plot of the condensate/superfluid density as a function of the tunneling parameter 
t at a fixed value of the chemical potential fi for each theory can be seen in Fig. 2. It 
is interesting to note that while the superfluid density from the field-theoretic approach 
increases linearly with t, the mean- field superfluid density quickly begins to fall off as t 
increases, a behavior which is at odds with the notion of a superfluid [13]. Furthermore, it 
can be seen that Eq. (102) is simply a first order series expansion of Eq. (103) about t = tc, 
meaning \ip\lq is the tangent line to I^M^f at t = tc- Thus, although the two results agree 
near the phase boundary, the mean-field prediction quickly begins to exhibit an unphysical 
behavior, suggesting that our field-theoretic result has a larger range of validity. 




(104) 




(105) 
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FIG. 3: Plot of the average number of particles per site as a function of the chemical potential fi/U 
with fixed t/U = 0.025. Left shows the mean-field prediction, while right shows the field-theoretic 
prediction. The solid blue curves are T/U = 0, and the dashed red curves are T/U = O.l/ks- 
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FIG. 4: Plot of the compressibility nU as a function of the chemical potential fx/U with fixed 
t/U = 0.025. Left shows the mean-field prediction, while right shows the field-theoretic prediction. 
The solid blue curves are T/U = 0, and the dashed red curves are T/U = O.l/ks- 



Next, we compare the average number of particles per lattice site (n) as computed in the 
two theories. In the ordered phase, the field theoretic prediction for (n) is given by Eq. (64) 
above, while the mean-field result is given by 

1 dTMF 



in) 



MF 



(106) 

IV'|'=IV'Imf 



In the Mott phase, = I'^Imf ~ 0, and both theories predict 



1 arp 

in) = Wmf = (107) 

which is simply {n)^. Plots of these two quantities as a function of the chemical potential 
at a fixed value of the hopping parameter are shown in Fig. 3. Although the predictions 
of both theories agree in the immediate vicinity of the phase boundary, we see that at 
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FIG. 5: Contours of constant (n) in parameter space for both the mean- field and effective action 
theories at T = 0. The dotted blue curve shows the first three lobes of the phase boundary. The 
solid blue curves are the predictions of the effective action theory, while the dotted red curves show 
the predictions of the mean-field theory. 

low temperatures {n)y^-p is not a monotonically increasing function of /i. In fact, the plot 
shows that away from the phase boundary, {n)y^^ actually decreases with increasing /i, i.e. 
the compressibility is predicted to be negative in the superfluid phase. This is directly 
at odds with the fact that a superfluid has, by definition, a positive compressibility [37]. 
The behavior of (n) as derived from the effective action, on the other hand, fits well with 
expectation further away from the phase boundary. This is highlighted in Fig. 4, which 
shows the compressibility k for each case. 

As a final point of comparison, we specify to the zero temperature case and examine 
contours in parameter space along which the average number density {n) is constant. Inside 
each Mott lobe, we know that the number density is fixed at the quantum number n of 
the lobe. In the superfluid phase, for fixed t we expect (n) to increase monotonically with 
increasing chemical potential. This implies that contours of constant {n) should be mono- 
tonic in /i. Such contours for (n) = 1, 2, 3 are shown in Fig. 5. Although the mean-field and 
effective action contours agree close to the lobe tip, the mean-field result exhibits a non- 
monotonic behavior farther away from the lobe tip. Given the non-monotonic behavior of 
(n)j^p discussed above, this is not wholly surprising. From these considerations, we conclude 
that the effective action theory has a larger range of validity than the mean field theory and, 
furthermore, we expect an increase in quantitative accuracy when higher powers of t are 
considered. 
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B. Effective Action Predictions in the Dynamic Case 



We now turn our attention to the Green's function and the zero-temperature excitation 
spectra of the Bose- Hubbard model. Using Eq. (77), we see that to first order in t the 
Green's function of the system in the ordered phase can be written as 

0'2 \^m) 



QiuJra, k) 



2af\a;„,0|0,cj^)(a™(0)) 



(0), 



6i- 



(4°>(a;„))4'" (0,0|0,0) 

where d is the lattice spacing, and k is restricted to the first Brillouin zone. At the phase 
boundary we have 6t — l/aj^^O) = 0, and the Green's function reduces to 



Q{uJm, k) 



(109) 



1 - 2ta'2\uJm) cos{k^d) 
This is precisely the result obtained in the Mott phase via a resummation of zero loop 
diagrams [32]. 

From Eq. (84), we see that the zero-temperature dispersion relation for amplitude exci- 
tations ujA{k) satisfies 



2aS(a;A,0|0,a;A)(4"i^(0)) 



(0), 







+ 



Qt + 



4*^(0) 



2t ^ cos (A;^d) , (110) 



4SK) ' (42(u;A))2aS(0,0|0,0) 
which we recognize also as the condition for poles in Eq. (108) continued to real time. 
While too complicated to solve exactly, Eq. (110) can be inverted numerically to yield the 
dispersion relation ujA{k). A plot of ujA^k) taken along the (1,1,1) direction in the first 
Brillouin zone is shown in Fig. 6. We observe that in the superfluid phase, WA{k) is gapped 
and quadratic, 

UAik) ^ A + r]P. (Ill) 

Furthermore, at the phase boundary, we find that the dispersion becomes gapless and linear. 

Next, we consider the zero-temperature dispersion relation u!g{k) of phase excitations. 
For the Bose-Hubbard model, Eq. (87) takes the form 

242(0)^ + 6t4°i(O) 







4°rM 4r(o) 42(0,010,0) 

+ b{0, 0, Ug, -Ug) - 2h{ug, 0, Ug, 0) - 2h{ujg, 0, 0, LOg 



26(0, 0, 0, 0) + 6(cu,, -cue, 0, 0) 
3 — cos {kfjd) 
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;il2) 



25 




FIG. 6: Plots of the zero-temperature dispersion relations tOAik) (Left) and uJe{k) (Right) for 
various values of the hopping t with fixed n = 1, /i/C/ = \/2 — 1 and with k = (1, 1, l)fc/\/3. The 
solid blue lines corresponds to t = tc ~ 0.028 [/, which for these values of /i and n is at the tip of 
the first Mott lobe. The dotted yellow lines corresponds to t = 0.03 U , and the dashed red lines 
corresponds to t = 0.035 U . Note that amplitude excitations exhibit a t-dependent energy gap, 
while the phase excitations are gapless in accordance with Goldstone's theorem. 



We can numerically solve this equation for ujg{k). A plot of uJo{k) along the (1, 1, 1) direction 
in the first Brillouin zone is shown in the right of Fig. 6. We see that in the superfluid phase 
the dispersion is quadratic with 

uJe{k)^Ck'^. (113) 

Finally, by comparing Eqs. (110) and (112), we find that at the phase boundary ua and cug 
are degenerate. Lastly, using the result (89), we can investigate the behavior of superfluid 
second sound excitations. From our observations above, we find a linear dispersion for small 
k, 

ujs(k) ^ ck (114) 

with the velocity of sound 

c=v^. (115) 

Thus, the velocity of second sound at any point near the phase boundary can be found 
from the above numerical inversions. With this, we plot the sound velocity c as a function 
of the tunneling t/U in Fig. 7. We observe that, since near the phase boundary a large 
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FIG. 7: Plot of the second sound velocity c as a function oit/U ioi t > tc ^ 0.028 U with n=l and 
li/U = V2-l. 

quadratic term suddenly appears in the phase dispersion relation, the sound velocity jumps 
immediately inside the superfiuid phase. This jump must be viewed cautiously, however, 
as the Ginzburg-Landau expansion is incapable of accurately describing critical behavior in 
the immediate vicinity of the phase boundary [27, 28, 41, 42]. As the mass of the phase 
excitations begins to increase faster than the gap in the amplitude excitations, we see that 
the sound velocity begins to decrease. Far from the phase boundary, however, we know 
from the seminal Bogoliubov theory that the sound velocity must increase as ^^t [36, 43], 
confirming that our theory is not valid in the deep superfiuid phase. 

VII. SUMMARY AND CONCLUSION 

In this paper we derived, to first order in the tunneling, the Ginzburg-Landau expansion 
of the effective action for a very general Bosonic lattice Hamiltonian. From the effective 
action we calculated many static and dynamic system properties of experimental interest. 
In specifying these results to the Bose-Hubbard model, we compared them with the corre- 
sponding findings of the standard mean-field theory. Although both approaches yield - up to 
first order in the tunneling - the same phase boundary, our method gives qualitatively bet- 
ter results deeper in the superfiuid phase. Additionally, we were able to find the dispersion 
relation for superfiuid excitations, which cannot readily be done in the mean-field approach. 
The primary advantage of our effective-action theory, however, lies in its extensibility. It is 
straightforward to generalize the derivation given in Sections II and III by calculating dia- 
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grams beyond the tree level in order to include higher-order tunneling corrections. As seen 
in Section V, this gives a systematic hopping expansion of the self-energy function in both 
the ordered and non- ordered phases, providing an arbitrarily precise description of the sys- 
tem dynamics near the phase boundary. This should allow, for instance, for the calculation 
of time-of-flight absorption pictures and their corresponding visibilities for the whole phase 
diagram [35]. Furthermore, given the generality of the formalism, our effective action the- 
ory can, in principle, incorporate a variety of interesting effects, such as disordered lattices 
[4, 26], vortex dynamics [36], and tunneling beyond nearest neighbor sites. In particular, 
an effective action for the disordered Bose-Hubbard model could give new insight into the 
nature of the Bose glass phase as a state of short-range order [13]. 
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